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Zooplankton (Cladocera, Copepoda) dynamics in the River Danube 
upstream and downstream of Budapest, Hungary 


Cs. VADADI-FULOP! 


Abstract. The spatial distribution and seasonal dynamics of zooplankton (Cladocera, Copepoda) were studied in the River 
Danube near Budapest, Hungary. The investigated river section was relatively poor in plankton, nauplii dominated. A total of 36 
species was recorded of which Acanthocyclops robustus, Thermocyclops crassus, Bosmina longirostris were the most abundant. 
There was a downstream increase in copepod densities, however, no other remarkable differences could be observed between the 
profiles upstream and downstream of the capital. Generally, the streamline was characterized by lower densities and lower 
number of taxa as compared to the river bank; nevertheless, there were differences between the left and the right banks both 
upstream or downstream as well. Seasonal dynamics was defined by a marked late winter-spring aspect and abundance peaks 
were found to be characteristic for Danube with high densities in May-June and August-September. Examining the relationship 
between zooplankton density and the hydrological regime, it can be concluded that zooplankton production in the main channel is 
of minor importance, rather floodplain areas and adjacent water bodies seem to be important sources of plankton biomass. 


pie Danube, as being the most species rich 
iver in the Palaearctic (Naidenow, 1998), is 
entitled to great interest, so the research of its 
plankton means a crucial challenge. The first in- 
vestigations on the zooplankton of the Hungarian 
section of the River Danube started in the 1910s 
with crustaceans (Kottász, 1913; Jungmayer, 
1914; Unger, 1916) and continued in the middle 
of the century (Woynárovich, 1944; Éber, 1955; 
Ponyi, 1962) mainly with descriptive works. Ne- 
vertheless, these studies served as a basis for 
further research. Dudich (1967) gave detailed in- 
formation about the fauna of the river Danube. 
Bothár has broken new ground in the studies of 
riverine plankton since the 1960s. With research 
among others on the sampling strategies (Bothár, 
1996), on the long-term status of planktonic 
crustaceans at the area of Göd (Bothár, 1968, 
1972, 1988a, 1994, 1996) and on the crustacean 
plankton near Budapest (Bothár, 1978, 1988b), 
she has made a major contribution to our better 
understanding of zooplankton in the River 
Danube. 


In the 1970s and particularly in the 1980s crus- 
tacean abundance increased notably, which was 
best professed by the abundance peaks in summer. 
The change in species composition and individual 
numbers implied eutrophication in Danube. How- 


ever, the hydrological regime has changed as 
well, which may influence zooplankton assem- 
blages both qualitatively and quantitatively (Bo- 
thár, 1985). Towards the estuary, crustacean den- 
sity increased by 25 fold within the 417 km long 
Hungarian river section (Bothár, 1988 a). The 
third period of the Danube plankton research be- 
gan in the 1990s with the studies of Gulyás (1994, 
1995, 1997, 2002) who has extended the research 
to rotifers. Rotifers dominated the plankton, only 
copepodites and nauplii were represented in si- 
milar abundance. 


The most frequent species were characteristic 
of eutrophicated stagnant waters and rivers of low 
current velocity (Gulyás, 1994, 1995, 2002). 
Large abundance of some species and the down- 
stream increase of zooplankton abundance im- 
plied increasing eutrophication in the river (Gu- 
lyás, 1995). Biomass values measured upstream 
and downstream of Budapest (near Budapest) 
were similar, then increased southwards (Gulyás, 
1997). 


Although numerous studies have been con- 
ducted on the zooplankton in the River Danube, 
relatively little emphasis was put on the transver- 
sal distribution of plankton (Bothár, 1978, 1985; 
Naidenow, 1971, 1979). On the other hand, de- 
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tailed surveys have been performed on the hori- 
zontal distribution of zooplankton along the river 
(Bothár, 1973, 1982, 1988 a; Gulyás, 2002; Pujin 
1990; Naidenow & Schewzowa, 1990, 1991). 


The objectives of this study were (1) to get a 
comprehensive picture of the crustacean plankton 
near Budapest, (2) to describe the community 
structure both from qualitative and quantitative 
points of view, (3) to present the spatial and tem- 
poral changes of microcrustacean plankton consi- 
dering the hydrological regime as a possible driv- 
ing force. 


MATERIALS AND METHODS 


Study sites 


River Danube is the second longest river in 
Europe, it is more than 2800 km long with a 
catchment area of 817,000 km?. The Hungarian 
section occupies 417 km. The middle section is 
377 km long and has a relatively balanced hyd- 
rological regime, which is controlled by the water 
discharge of the upper section. In the middle 
section, current velocity ranges between 0.8—1.2 
m sec” in the upper layers, which can reach the 
values of 2-2.5 m sec’ during flooding events. 
The water level fluctuation takes 5-8 m, discharge 
values range from 500 (low water period) to 
6000-8000 (flood) m? s”. The shoreline is mostly 
regulated, rip-rap is characteristic of this section. 


The present study was conducted in the main 
channel of the River Danube at two profiles up- 
stream and downstream of Budapest (Üjpest—Bé- 
kásmegyer, 1657 rkm and Tókól-Százhalombatta, 
1623 rkm) (Fig. 1). The general characteristics of 
the two profiles are similar, but the sampling site 
of Tókól-Százhalombatta is characterized by gra- 
vel banks interrupted occasionally with rip-rap, 
whereas in the profile of Üjpest-Békásmegyer 
rip-rap is the main component. There was no mac- 
rovegetation either at the upstrem or at the down- 
stream profile throughout the years. The sampling 
point of Békásmegyer is situated in the down- 
stream section of the Danube arm of Szentendre, 
near to the estuary into the main channel. 
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The study sites were designated on the basis of 
our objectives, i.e. possibilities for transversal 
sampling (ferry), and sampling profiles upstream 
and downstream of Budapest. 


Sampling and data analysis 


Samples were collected at biweekly intervals 
from October 2006 to September 2008 at two 
sampling profiles each containing three sampling 
sites (stream line, left bank and right bank). 
Samples taken at the river bank were performed at 
water depth of some 1-2 metres. During the 
winter period (between December and February) 
zooplankton was sampled monthly. 100 litres of 
water were taken from the surface water layer 
(upper 50 cm water column) and filtered through 
a plankton net (50 um mesh size). The material 
collected was preserved in situ in 496 formal- 
dehyde solution. A total of 270 samples were col- 
lected and analysed. Nauplii were counted in 5 ml 
subsamples in special counting chambers after 
homogenization. For the taxonomic determination 
of the animals identification keys by Gulyás and 
Forró (1999, 2001), Einsle (1993), Amoros (1984) 
and Dussart (1969) were used. Copepods and cla- 
docerans were identified to species level, how- 
ever, copepods belonging to the suborder Har- 
pacticoida and ostracods were only counted. 


Water temperature was measured in situ 
whereas conductivity was measured in laboratory. 
Water discharge and water level data (1645.5 
rkm) were obtained from the Environmental and 
Water Research Institute (, VITUKT"). Water re- 
sidence time was calculated with the formula: 


R = 0.08A4^/Q"' 


where R is the residence time at the sampling site 
(days), Aq watershed area upstream of the sam- 
pling site (km?), and Q river discharge (m? s”) 
(Soballe & Kimmel, 1987). We used the dis- 
charge data measured at Budapest (1646.5 rkm). 
In order to explore the temporal patterns, cluster 
analysis and non-metric multidimensional scaling 
(NMDS) using the Euclidean distance were per- 
formed (with standardized data). 
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Correspondence analysis was applied to de- RESULTS 


monstrate the relationship between species and 
sampling sites, that is certain species seem to be 
more abundant at a definite sampling site. Di- 
versity t-test 1s a useful tool to compare the di- 
versities of the sampling sites by calculating the 
Shannon diversity of each sampling site and com- 
paring the diversities statistically. Linear corre- 
lation was used to detect any significant asso- 
ciation between some environmental variables 
and zooplankton community. When necessary, 
data were transformed (log, square root) to obtain 
the normal distribution. Significant differences 
were identified at p « 0.05. All data analyses were 
performed using the PAST program (Hammer ef 
al., 2001). 


During the study period, water level varied 
between 83 and 689 cm-s, maxima were recorded 
in September of 2007. This marked peak was two 
times higher than other peaks. Low water level 
occurred in the year of 2006, whereas relatively 
high water level persisted in spring and summer 
of 2008. Discharge values (ranged between 1030 
and 6420 m? sec”) strongly followed this pattern. 
Despite some deviations, the annual discharge 
pattern typical of river systems in the temperate 
region was observed with low water period in 
autumn and winter, and discharge peaks in early 
and late spring. Water residence time ranged 
between 48 and 58 days (Figs. 2 a-d). Water 
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Figures 2 a-c. Hydrological parameters recorded during the study period. a 


90 


Vadadi-Fülóp: Zooplankton dynamics in the River Danube 


Water temperature (^C) 


10 


2006.10.03 
2006.11.03 
2006.12.03 
2007.01.03 
2007.02.03 4 
2007.03.03 
2007.04.03 
2007.05.03 4 
2007 06.03 
2007.07.03 
2007.08.03 4 


2007.09.03 


2007.10.03 


2007.11.03 4 
2007.12.03 
2008.01.03 
2008.02.03 4 
2008.03.03 
2008.04.03 
2008.05.03 
2006.06.03 
2008.07.03 
2008.08.03 4 
2008.09.03 


Figure 2 d. Water temperature recorded during the study period 


temperatures showed typical annual cycle, ma- 
ximum values were observed in June-August 
(23-24?C) and minimum in January (0.5 °C). 
There was a negative relationship between water 
level and conductivity both upstream (r = 0.40; p 
= 0.0073) and downstream (r = 0.38; p = 0.0118) 
of Budapest over the study period. 


A total of 36 species were observed of which 
11 belonged to copepods and 25 to cladocerans 
(Table 1). Diacyclops crassicaudis was recorded 
for the first time in the Hungarian section of the 
River Danube. The occurrence of Diaphanosoma 
mongolianum in the main channel of the Danube 
is new, as it was only reported from the Szigetkóz 
backwaters and few water bodies respectively 
(Gulyás & Forró, 1999). Frequent species includ- 
ed Acanthocyclops robustus, Thermocyclops cras- 
sus, Bosmina longirostris, Alona rectangula, Chy- 
dorus sphaericus, moreover harpacticoid cope- 
pods were also determining. The species number 
in the streamline was similar upstream and down- 
stream (17-16 species), 10 common species were 
recorded. The distribution of the species number 
within the sampling profile was more dispersed 
downstream with a maximum of 30 species de- 
tected at Tököl. The relative abundance of nauplii 
was mostly greater than that of other zooplankton 
groups, whereas adult copepods and copepodites 
were represented in similar abundances. Clado- 


cerans contributed up to 19.596 of the total density 
upstream of Budapest and 18.1% downstream, 
respectively. 


The percentage compositions of the zooplank- 
ton community upstrem and downstream of the 
capital are presented in Figs. 3—4. There is evi- 
dence of a remarkable downstream increase in 
copepod densities meaning that adult copepod 
densities became doubled and copepodit densities 
increased by 18%. This phenomenon is mainly 
due to the fact that the relative contribution of 
Thermocyclops crassus increased considerably at 
the downstream profile, although calanoid den- 
sities increased as well. Such notable longitudinal 
differences cannot be seen either in cladoceran 
community composition, or in nauplii. 


Generally, the streamline was characterized by 
lower individual numbers and lower number of 
taxa in comparison with the river bank (Table 2), 
however, there were also differences between the 
left and the right banks both upstream or down- 
stream of Budapest. Average densities ranged be- 
tween 5.08 and 12.96 ind./100 litres, whereas 
maximum density reached 87 ind./100 litres at 
Békásmegyer. Note that nauplii densities are in- 
cluded as well, these results suggest that the river 
is poor in planktonic crustaceans. Regarding di- 
versities it seemed to be essential to present the 
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Table 1. Zooplankton taxa recorded in the River Danube (October 2006 to September 2008) and their codes. Abbreviations: B: 
Békásmegyer (right bank); UBS: Ujpest-Békásmegyer streamline; U: Ujpest (left bank); Sz: Százhalombatta (right bank); 
TSzS: Tókól-Százhalombatta streamline; T: Tököl (left bank) 


Taxa Code! B UBS U Sz | TSzS T 
06/07/08/06/07/08/06/07]08/06/07]08/06/07]08/06/07/08 

CLADOCERA | | | a 
Alona affinis (Leydig, 1860) Aaff + + + + 
Alona guttata Sars, 1862 Agut + + 
Alona intermedia Sars, 1862 Aint + 
Alona quadrangularis (O. F. Müller, 1785) Aqua + 
Alona rectangula Sars, 1862 Arec | 
Bosmina coregoni Baird, 1857 Bcor + 
Bosmina longirostris (O. F. Müller, 1785) Blon 
Chydorus sphaericus (O. F. Müller, 1776) Csph | 
Daphnia cucullata Sars, 1862 Dcuc + H 
Daphnia longispina O. F. Müller, 1785 Dlon + F + 
Diaphanosoma brachyurum (Liévin, 1848) Dbra + + 
Diaphanosoma mongolianum Uéno, 1938 Dmon + F 
Disparalona rostrata (Koch, 1841) Dros H + 
Graptoleberis testudinaria (Fischer, 1848) Gtes + 
Iliocryptus sordidus (Lievin, 1848) Isor t| | [| H it 
Leydigia acanthocercoides (Fischer, 1854) Laca | | | [t| jJ 
Leydigia leydigi (Schoedler, 1863) Lley + + + 
Macrothrix hirsuticornis Norman & Brady, 1867 |Mhir + H + is 
Macrothrix laticornis (Fischer, 1848) Mlat * + 
Moina macrocopa (Straus, 1820) Mmac| |+ + 
Moina micrura Kurz, 1874 Mmic | 
Pleuroxus aduncus (Jurine, 1820) Padu + + 
Scapholeberis mucronata (O. F. Müller, 1785) Smuc| |+ BS + + I || 
Sida crystallina (O. F. Müller, 1776) Scry || + | 
Simocephalus vetulus (O. F. Müller, 1776) Svet + 
COPEPODA 
CALANOIDA 
Eudiaptomus gracilis (Sars, 1863) Egrac + + 
Eurytemora velox (Lilljeborg, 1853) Evel + 
CYCLOPOIDA 
Acanthocyclops robustus (Sars, 1863) Arob + | 4 T H H 
Acanthocyclops vernalis (Fischer, 1853) Aver ||] | + | 
Cyclops strenuus Fischer, 1851 Cstr EN ++! 
Cyclops vicinus Uljanin, 1875 Cvic + 
Diacyclops crassicaudis (Sars 1863) Dcra * 
Eucyclops serrulatus (Fischer, 1851) Eser | H 
Mesocyclops leuckarti (Claus, 1857) Mleu + + 
Paracyclops fimbriatus (Fischer, 1853) Pfim + + 
Thermocyclops crassus (Fischer, 1853) Tcra 
HARPACTICOIDA Har |+ | 
OSTRACODA Ost + 
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Table 2. Major ecological parameters of the zooplankton community (n — 44). Density data include nauplii as well. For the 
abbreviations of the sampling sites see Table 1 


Sampling site B UBS U Sz TSzS T 

Average density (ind./100L) 8,48 6,80 9,82 9,17 5,08 12,96 
Maximum density (ind./100L) 87,00 49,00 43,00 46,00 22,00 60,00 
Taxa S 26,00 17,00 25,00 19,00 16,00 30,00 
Dominance D 0,23 0,13 0,08 0,11 0,13 0,14 
Shannon H 2,21 2,31 2,83 2,52 2,32 2,48 
Evenness e^H/S 0,35 0,59 0,68 0,66 0,64 0,40 
Berger-Parker 0,45 0,26 0,16 0,19 0,21 0,25 
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Figure 3. The percentage composition of the zooplankton community upstream Budapest. Results include data obtained from the 
streamline, left and right river banks (summed) 


results of two different diversity indices. Shannon 
diversity values ranged within a relatively close 
interval (2.21-2.83), highest diversity and eveness 
values were measured at Ujpest. Berger-Parker 
diversity peaked at Békásmegyer where the Shan- 
non diversity was the lowest, which is due to the 
dominance of some species (dominance D - 


0.23). Based on the diversity t-test, sampling 
station Ujpest (U) differed significantly from the 
others (U-B: t = 3.89, p<0.001, U-UBS: t = 3.77, 
p<0.001, U-Sz: t= 2.29, p = 0.02318, U-TSzS: t = 
3.68, p<0.001, U-T: t = 2.54, p = 0.0156). One 
more significant difference included Bekäsme- 
gyer-Szäzhalombatta (t = 2.32, p = 0.0214). 
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The correspondence analysis (Fig. 5) suggests 
that the sampling site Békásmegyer, which is 
situated in the side arm of Szentendre, has spe- 
cific features. The high individual numbers of 
Moina micrura and the presence of Simocephalus 
vetulus, Alona guttata, Mesocyclops leuckarti, 
Graptoleberis testudinaria are responsible for its 
separation. The sampling site Ujpest (U) is also 
determined by the presence of several taxa. The 
sampling profile downstream of Budapest seems 
to form one group, although the sampling site 
Tököl (T) has some characteristic elements, e. g. 
Bosmina coregoni, Sida crystallina and it is do- 
minated by Thermocyclops crassus. 


Examining the temporal patterns of the zoo- 
plankton assemblage, no clear trajectory could be 
detected, however, there was a marked late win- 
ter-spring aspect (marked) with great similarities 
among samples. Samples taken in summer and 
autumn form a second group, but are far more 
dispersed (Fig. 6). Species number was low dur- 
ing the winter and increased in late spring, then 
remained relatively high during the summer 
months and decreased in autumn. Total abun- 
dance, however, peaked in spring of 2008 and in 
summer of 2007. The seasonal dynamics can be 
characterized by two population peaks within the 
year, the first in May-June, the second in August- 
September. 


The seasonal patterns were similar between the 
upstream and downstream sections, with a posi- 
tive correlation of the zooplankton community 
density between the two sampling profiles over 
the study period (r = 0.58; p«0.001). A strong po- 
sitive correlation was detected between water 
temperature and zooplankton density at Tököl- 
Százhalombatta (r = 0.35; p = 0.0178), but not at 
Üjpest-Békásmegyer, however, zooplankton den- 
sity without nauplii was positively related to wa- 
ter temperature (r = 0.42; p = 0.004). Zooplankton 
density (without nauplii) was positively correlated 
with water level (r = 0.32; p = 0.0303) and 
discharge (r = 0.32; p = 0.0307) and negatively 
correlated with water residence time (r = 0.32; p = 
0.0305) at the upstream section. However, these 
relationships were not significant at the down- 
stream profile. 
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DISCUSSION 


The time available for zooplankton to develop 
was esimated at 48—58 days (an estimate of the 
time the water has been in the river system) 
assuming that plankton drifts with the current 
passively. During this time 5 generations of cla- 
docerans and 1 or 2 generations of copepods may 
develop (Naidenow, 1998). The growth is sup- 
posed to depend on the hydrological regime (in- 
creased residence time favours zooplankton, that 
is zooplankton benefits from low water velocity). 
This 1s not only due to the mechanical effect of 
the drift, but it is connected with the fact, that 
suspended matters have a negative impact on 
planktonic crustaceans (Zsuga et al., 2004; Gu- 
lyás, 2002;). However, we have found a negative 
relationship between water residence time and 
zooplankton density and a positive relationship 
between water discharge and zooplankton density, 
which strongly suggests that zooplankton pro- 
duction in the main channel is of minor impor- 
tance, which is consistent with the findings of 
Reckendorfer et al. (1999). Floodplain areas and 
adjacent water bodies seem to be rather important 
sources of plankton biomass (Saunders & Lewis, 
1989; Naidenow 1998; Schiemer et al. 2001; Zsu- 
ga et al., 2004). The positive correlation of zoo- 
plankton density with water level also seems to 
support this hypothesis. According to Recken- 
dorfer et al. (1999), the physical interaction of 
flow regime and river margin morphology de- 
termines the availabihty of inshore storage zones 
and the rate at which plankton are added to the 
main river channel. 


The dominant species that were found in the 
present study in the main channel (Acanthocyc- 
lops robustus, Bosmina longirostris) are consis- 
tent with the findings of Bothár (1985, 1988 b, 
1994) and Gulyás (1994, 1995, 2002), however, 
the relatively large contributions of Thermocyc- 
lops crassus and Harpacticoida are new. Two 
thermophyl species, Thermocyclops crassus and 
Moina micrura are reported to become abundant 
in the River Danube since 1971 (Bothár, 1975), 
but they were regarded only as secondary species 
at Gód (rkm 1669) (Bothár, 1985). The relatively 
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Figure 4. The percentage composition of the zooplankton community downstream Budapest. Results include data obtained 
from the streamline, left and right river banks (summed) 


large frequency and abundance of Thermocyclops 
crassus may be related to the changing hydro- 
logical regime and increased temperature. Beyond 
the above-mentioned, there is no evidence of any 
significant long-term changes in species com- 
position. The population peaks observed in the 
present study (May-June, August-September) are 
in line with the results of Bothár (1985, 1978), 
Bothár and Kiss (1990) and Gulyás (1995). Re- 
garding densities, the examined river stretch was 
poor in plankton, the individual numbers were 
similar to the findings of Bothár et al. (1971) and 
Bothár (1978) at Újpest-Békásmegyer in the 
1970s, at Adony (1598 rkm) in the 1980s (Bothár, 
1988b), but not to that of Göd (Bothár, 1994). In 
quantitative aspect, contrasting my results to that 
of Gulyás (1994, 1995, 1997) is not meaningful 
since the author presented the overall densities 
(with rotifers). 


There were no remarkable differences between 
the upstream and downstream sections, only co- 
pepod density increased downstream. This does 


not support the hypothesis that municipal and in- 
dustrial wastewaters of the capital may create sig- 
nificant impact on zooplankton assemblages 
downstream Budapest. Bothár (1988) has observ- 
ed lower maximum densities and different sea- 
sonal dynamics patterns downstream of Budapest. 
According to the author, these differences can be 
explamed by the effects of wastewaters and 
changed stream conditions due to the regulated 
river bed and the side arm of Szentendre. Ne- 
vertheless, Gulyás (1997) has not found any 
remarkable differences between the upstream and 
downstream sections. Although the examined ri- 
ver stretch is poor in planktonic crustaceans, the 
recorded species and their contributions to the 
total density suggest that the river has a moderate- 
ly rich fauna with several tichoplanktonic ele- 
ments. The latter also implies for the significance 
of inshore storage zones. 


The spatial distribution of zooplankton across 


the river was not equal, the streamline was cha- 
racterized with lower densities and lower number 
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Figure 5. Biplot of the correspondence analysis. Codes of the taxa are presented in Table 1. Abbreviations of the sampling sites: 
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Figure 6. NMDS plot of the samples. Numbers represent the months: 1-January; 2-February etc. 
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of taxa. Similarly, densities were often found to 
be higher nearshore (Bothár, 1978, 1985; Thorp et 
al., 1994; Mitsuka & Henry, 2002). 


In summary, my results pointed out the evi- 
dence of spatial heterogeneity across the river, 
however, save adult copepods no major difference 
was detected between the upstream and down- 
stream sampling profiles. The author stresses the 
importance of adjacent lentic areas as sources of 
planktonic crustaceans. 
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